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^ The  transient  response  of  s system  of  two  concentric  spherical  elastic  shells  coupled  by  an 
ideal  fluid  and  imping^  by  an  incident  plan  pressure  pulse  is  analyzed.  The  classical  techniques 
of  separation  of  variables  and  Laplace  transforms  are  employed  for  solving  the  wave  equations 
governing  the  fluid  motions  and  the  shell  equations  of  motion.  A scheme  of  iterative  convolution 
was  devised  for  the  inversion  of  the  Laplace  transforms  that  facilitates  the  calculation  of  accurate 
transient  solutions  of  the  response  of  the  shells.  A sample  calculation  of  shell  responses  was 
performed  and  results  are  compared  to  the  case  In  which  the  outer  shell  is  absent.  This  set  of 
results  demonstrates  that  a thin  outer  shell  tends  to  be  transparent  to  the  incident  pulse.  ^ 
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TRANSIENT  RESPONSE  OF 
TWO  FLUID-COUPLED  SPHERICAL  ELASTIC  SHELLS 
TO  AN  INCIDENT  PRESSURE  PULSE 


INTRODUCTION 

An  underwater  weak  shock  wave  sufficiently  far  away  from  its  generating  explosion 
source  is  often  treated  as  an  acoustic  pulse.  There  is  prolific  literature  on  the  studies  of  the 
transient  interaction  among  such  pressure  waves  and  single  elastic  shells  of  simple  shapes  (IJ). 
The  results  not  only  reveal  many  of  the  essential  physical  phenomena  involved  in  the  interac- 
tion problem  but  also  are  quite  useful  for  the  verification  of  approximation  methods  for 
predicting  the  underwater  explosion  response  of  submerged  structures  surrounded  by  an  exte- 
rior fluid  medium  of  infinite  extent  (2-51.  In  the  present  endeavor,  the  transient  response  of  a 
system  of  two  fluid  coupled  concentric  spherical  elastic  shells  impinged  by  an  external  incident 
plane  shock  wave  is  analyzed.  This  purports  to  gain  physical  insight  in  the  response  as  well  as 
to  provide  a data  base  for  the  development  of  general  numerical  methods  for  predicting  the 
underwater  explosion  response  of  fluid  coupled  shell  systems  such  as  the  double  hull  section  of 
a submarine.  The  problem  with  the  spherical  geometry  permits  the  separation  of  variables  in 
the  wave  equation  governing  the  fluid  motion  and  the  shell  equations  of  motion.  The  use  of 
Laplace  transforms  then  facilitates  the  calculation  of  satisfactory  transient  solutions  of  the 
response  of  the  shells. 

A paper  appears  in  the  Russian  literature  (6)  independently  dealing  with  exactly  the  same 
problem.  The  results  obtained,  however,  only  pertain  to  the  point-symmetric  and  the  transla- 
tional motions  of  the  shells  and  these  two  terms  of  the  series  solution  are  insufficient  for  the 
description  of  the  complete  shell  response.  It  would  also  seem  that  the  scheme  of  dual  Volterra 
integral  equations  for  the  calculation  of  the  inverse  Laplace  transform  used  in  Ref.  6 is  un- 
necessarily complex  and  numerically  inefficient.  It  can  be  readily  shown  that  only  a single  and 
simpler  integral  equation  is  needed  if  this  scheme  is  used. 

In  this  report,  a simple,  straightforward  method  is  devised  for  the  inverse  Laplace 
transform,  and  the  asymptotic  behaviors  of  the  shells  are  analytically  discussed.  An  example 
calculation  is  also  carried  out  using  eight  terms  of  the  series  solution  for  adequate  convergence 
of  shell  displacements,  velocities,  and  stresses.  Time  histories  of  the  transient  reponse  of  the 
interior  shell  are  presented.  This  example  demonstrates  that  a thin  exterior  shell  tends  to  be 
transparent  to  the  incident  pulse. 


ManuKript  submitled  December  16,  1977. 
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Description  of  the  Probiem 

Figure  1 sketches  the  fluid-coupled  spherical  shell  system  and  the  incident  plane  pressure 
wave.  The  fluid  surrounding  the  outer  shell  and  that  between  the  two  shells  are  considered  to 
be  ideal  compressible  fluids  in  linear  wave  motions  and  can  be  characterized  by  their  unper- 
turbed mass  densities  and  sound  speeds,  i.e.,  by  (p'.  c')  and  (p.  c),  respectively.  The  shells  are 
initially  concentric.  In  this  study,  the  strength  of  the  incident  wave  is  sufficiently  weak  such 
that  the  shell  deflections  are  elastic  and  small  and  the  deviation  from  the  concentricity  remains 
negligible  for  the  time  duration  of  interest.  The  mass  densities.  Young’s  moduli  and  Poisson’s 
ratios  of  the  outer  and  inner  shells  are  (p',  E'.v')  and  (p,.  £,  »>),  respectively.  The  middle 
surface  radii  and  thicknesses  of  the  outer  and  inner  shells  are  (o',  h')  and  (a,  A),  respectively. 
The  <t>-coordinate  of  the  spherical  system  (r,  0,  <t>)  is  not  shown  in  Fig.  1 since  it  is  not  needed 
due  to  the  symmetry  of  the  problem.  The  origin  O coincides  with  the  unperturbed  center  of 
the  shells. 
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Fit.  1 ~ Gcomcdy  of  the  problem 


The  deflections  of  the  inner  shell  in  the  r-  and  0-direction,  normalized  with  respect  to  the 
outer  shell  radius  o',  are  denoted  by  w and  u respectively  and  those  of  the  outer  shell  by  h' 
and  u'  respectively.  The  total  pressure  field  exterior  to  the  outer  shell  is  denoted  by  p'  (r,  0,  i) 
and  that  between  the  shells  by  p(r,  0,  t)  where  / designates  time.  The  following  dimensionless 
parameters  will  be  used  in  the  mathematical  formulation: 
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II ‘'and  II  satisfy  the  wave  equations 

V^ir  (2) 

and 

V^II-r;-^  (3) 

respectively,  where  V ^ is  the  Laplacian  operator.  The  boundary  conditions  of  the  problem  are 
that  11 ' satisfies  the  radiation  condition  at  far  field  and  that 

dll  _ dll' d^  ..  o _ . 


dR  dR 


R - I 


M -t.  (5, 

All  quantities  except  the  incident  pressure  field  have  quiescent  initial  conditions. 

A Laplace  transform  pair  is  defined  as 

m 

w(9.s)  - ^ w{».  T)e  -^dT 

y +>— 

w(9,T)  — — — f w{9,s)e'^ds  (6) 

27nJ 

where  y lies  to  the  right  of  all  singularities  of  win  the  complex  s-plane  and  i — ( — I ) 

Due  to  the  spherical  geometry  of  the  problem,  the  solutions  can  be  expanded  in  terms  of 
series  of  Legendre  polynomials  as  the  following. 

UiR.e.T)  - £ n„{R.T)P„{cos9) 

m H) 
oo 

ii'iR.e.  r)  - £ n'{/f,  T)p„kos9) 

m M) 

w(9,T)  - £ w^(T)P„{cos  9)  (7) 


w'ffl.  n - £ w^{T)P„  {cos  9) 

m H) 

dP„  (cos  9 ) 

- I 

m -1 

~ dP„  (cos  9 ) 

— 

m-l 

where  P„  is  the  Legendre  polynominal  of  the  first  kind  and  mth  degree. 
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In  the  Laplace  transform  domain,  the  equations  of  motion  of  the  elastic  shells  are  17) 

-A/{IUU.5)({^5^  -L  g^) 


-A/Mll'd.s)  -ll,„(l.s)|(i^  +«') 

ykt“  M ■ I ■ ■ ■ " 

^ 4 • V J I r 

s + A's^  + 

m “ 0.  I,  2, ....  (8) 

({^S^  ■'"r+7“'" 

(s^  + a')?'  - cl  + - ^ K 

m - 1.2.3 (9) 

It  should  be  noted  that  the  solution  method  developed  here  is  applicable  for  any  linear  elastic 
theory  for  the  spherical  shells.  The  choice  of  the  version  in  Eqs.  (8)  and  (9)  is  for  the  compari- 
son of  results  previously  obtained  in  Ref.  7. 

The  total  pressure  field  exterior  to  the  outer  shell  II 'consists  of  the  pressure  due  to  the 
incident  wave  and  those  due  to  scattering  and  radiation  by  the  outer  shell  An  arbitrary  in- 
cident plane  pressure  wave  impinging  the  vertex  of  the  outer  shell  (R  * l.w  “ 0)  at  T ™0 
can  be  expressed  by  the  following  series  [7] 

ee 

ir{/?.».5)  £ (2m  -I-  l)/„,(/li)/»„(cosfl).  (10) 

m ^0 

where  /(s)  is  the  Laplace  transform  of  the  time  characteristics  of  the  incident  wave  and 
i„(Rs)  is  the  abbreviated  notation  for  the  modified  spherical  Bessel  function  of  the  first  kind 
[W(2Rs)1'/2/^  + ,/2(/?s)  (8). 

Solutions  In  the  Laplace  Transform  Domain 

It  can  be  shown  that  the  solutions  to  the  system  of  Eqs.  (2)  through  (10)  are: 

lT'(/?.s)  - (2m  -t-l)  - [i„(Rs)k„{s)  - k„(Rs)i'„(s)] 

k„  (s) 

5w'X„(/?5) 

km(s) 


n„(/?.  s) ^ ^ ^ : 

fr  I'm  “ 'm  (frC*)  1 
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-^/(i)A/^(s^  +a;)|((V  +X«(V  +M«) 

lC(v)C(c,{i)  -C(c,{s)*„(c,s)l  +P,A/-^  ((^s^  +a„) 

X l/„(c,5)*„(c,{s)  - /„(c,(s)*^(c,i)l)  (13) 

A * PrMM'U  +o„)(s^  +a').  (14) 

2s^  2c/{5 

In  the  above  equations  k„{s)  is  the  abbreviated  notation  for  the  modified  spherical  Bessel 
function  of  the  third  kind.  (ir/2jl  + i/j  (s).  The  prime  denotes  diflerentiaiion  of  the 
Bessel  functions  with  respect  to  their  arguments,  and 

+ Mm)l*;i,(c,s)C(c,{i)  - k„(c,is)i„(c,s)\ 

+ a„)I/^(f,i)A„(f,(5)  - /f^(c,5)»^(c,(s)l| 

X Ik;,(s)(s^  + X;;s2  +m;)  - W'sX„(s)(s2  + a')l 

+a^)(i*s*  + K„ih^  +p„)k'„(s) 

X (x„(c,{s)/„(c,j)  - /„(c,{s)X,„(c,s)I 

Ip  |2 

+ -^  MM'ishs^  + + a„)k„{s) 

X I/„(c,j)X„(c,{s)  - /„(f,{j)*„(f,s)l  (15) 

With  us:  of  the  Tauber's  theorem  of  Laplace  transforms  (9],  some  of  the  asymptotic 
behaviors  of  the  shell  responses  at  late  time  can  readily  be  revealed  from  Eqs.  (13)  and  (14). 
Specifically,  for  the  case  where  the  incident  wave  n ' is  a unit  step  wave.  i.e../(s)  - 1/s, 

UT)  -2CV(1  +ii^p, MM' 

4C^C^c}{\  -{^)  +6p,Mi^C}  6p,M'C^ 


wo(T) 


r-»  4C^C,^c/(l  -{^)  +6p,Mi^C^  + 6p,M'C^' 

These  shell  deflections  occur  long  after  the  incident  wave  has  engulfed  the  outer  shell  and  can 
also  be  found  by  static  analysis. 


Again,  for  the  unit  step  incidence  case, 


wf  (T) 

f-*oo 


-3A/'((1 +2{^)p,A/+6(l  -{^)1 


6p,A/(l  +2C^)+6yl/M(l  -C^)  + (2 +{^)p,J +P,A/A/'(1  +2(^+2p,(l  -(^)]+36(l  -{^) 
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and 

wi  (T)  - 
r-- 

-9p,MM' 

6p,A/(l  +2{^)+6A/M(l  -{^)  + (2+{^)p,l+p,A/A/'Il  +2{^+2p,(I  -{»)l+36(l  -{M 

(19) 


where  the  dot  denotes  differentiation  with  respect  to  T.  These  are  late  time  translational  veloci* 
ties  of  the  shells  in  the  direction  of  propagation  of  the  incident  wave.  It  can  be  seen  that  the 
only  condition  under  which  Wj  -•  at  late  time  is  p,M  ~ 3,  i.e.,  when  the  inner  shell  is 
neutrally  buoyant  in  the  interior  fluid.  Otherwise,  they  are  not  equal  to  each  other.  Therefore, 
when  interpreting  the  results  of  the  present  problem  for  large  values  of  T,  the  relative  posi- 
tions of  the  shell  must  also  be  examined.  If,  at  some  time,  the  deviation  from  concentricity  be- 
comes excessive,  the  results  of  analysis  thereafter  are  no  longer  physically  meaningful. 

Formulae  (16)  through  (19)  are  convenient  for  parametric  studies  of  the  effects  of  vari- 
ous shell  properties  and  fluid  properties  on  the  symmetric  and  translational  responses  of  the 
shells,  and  they  are  also  quite  useful  for  providing  asymptotic  checks  for  the  numerical  calcula- 
tions. 

The  Inversion  of  the  Transformed  Solutions 


The  spherical  Bessel  functions  can  be  expressed  in  finite  series  of  elementary  functions 
[8].  If  this  property  is  used,  the  transformed  solutions,  e.g.,  Eq.  (14),  can  be  rearranged  in  the 
following  form: 


w„(s) 


(20) 


where 


r^(j)  - -2(-l)'"  + >(2/n  + 1)(pA,)/(s)A/A/'(C5)'"  + '(c,j)2'”+2(52  +a')(C^s^  + a„) 
“ lx„(s)x„(  -c,s)(s*  +a's^  +p')  + + a' ) I)'^  (s)3r„  ( -c,j) 

+ p,MCh^(Qh^  +ajr„(c,Cs)] 

- U„(s)-r„(c,s)(s^  +A's2  + + a;;,)lY„(s)X„(c,s) 

-p,l'^(c,s)jr^(s)l)(jr„(-c,(s)(£*5<  +X^£^j2 

+ p,Mi^sHih^  +a^)K^(-c,£5)l 

m + l m 

- L ('W  +3/2,k)2  -*5'"  + ' + 1/2,  k)2  -‘s'"-* 

*-o  * -0 


+ 1/2,  k)2  -‘s'"-* 

k -0 

(m  -I-  kV. 


(m  -I-  1/2,  k)  - 


k!(m  - A)!  • 


(21) 
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Both 
case,  r 


A„  and  B„  are  (3m  + ll)rA  order  real  polynomials  of  s and,  for  the  plane  incidence 


is  a (3m  + 6)/Aorder  real  polynomial. 


Equations  of  the  form  of  Eq.  (20)  appear  in  many  applied  problems  such  as  wave  propa- 
gation in  layered  media,  transmission  lines,  etc.  The  exponential  factors  in  the  denominator 
signify  the  reflections  of  waves,  in  the  present  case,  to  and  from  ( and  R * ! The  di- 
mensionless time  for  a disturbance  to  traverse  the  distance  (1  - ()  is  ij.  There  is  a variety  of 
schemes  for  calculating  the  inverse  Laplace  transform  of  Eq.  (20).  One  of  the  standardized 
techniques  is  the  D’Alembert  expansion  (10)  by  which  Eq.  (20)  can  be  rewritten  as 


W„(5)  - 


A„(s) 


1 + 


-2i?J 


-2<n|i  , 


-4»|i 


(22) 


n - 1,2,3.... 

Since  A„  is  a real  polynomial  of  s,  its  roots  can  be  accurately  computed  by  well-established  nu- 
merical procedures  and  the  inversion  of  every  individual  term  on  the  right-hand  side  of  Eq. 
(22)  can  then  be  obtained  by  the  method  of  residue.  The  calculation  of  residues,  however, 
would  be  rather  clumsy  for  terms  with  A„  of  high  power.  The  following  scheme  utilizing  the 
convolution  theorem  is  devised  to  circumvent  this.  From  Eq.  (22), 

**'„(?')  - H-°(r  -r,)H{T  -T,)  -I-  wl{T  -3t,)  + wl{T  - Sr,) H(T  - 5rj) 


-»■  •••  + w^lT  - On  + l)T|l//lr  - (2n  + l)i)l  + 
where  H is  the  Heaviside  step  function  and 

-|G„(r  -r)H-O(r)*. 

T 

-wim  - -r)wHr)dr. 


(23) 


n —I 


(n 


] 


— t)»v"  ’ (t)</t 


(24) 


In  Eq.  (24), 


and  G-(T)  are  the  inverse  transforms  of  r_(s)/,4_(s)  and 


[A„  (s)  — B„  (s)l//4„  (s)  respectively.  They  are  to  be  first  accurately  obtained  by  the  method 
of  residue.  This  is  quite  practical  with  the  current  computing  technology  for  up  to  moderate 
values  of  m.  For  large  m asymptotic  expansions  for  the  spherical  Bessel  functions  can  be  used 
and  A„  and  B^  will  assume  simpler  forms.  For  the  calculations  of  shell  responses  in  the 
present  problem,  large  m terms  are  not  needed.  The  successive  convolutions  required  in  Eq. 
(24)  can  be  conveniently  programmed  and  carried  out  in  a modern  digital  computer.  Since 
w®  (D  and  G„(T)  are  composed  of  terms  formed  by  an  exponential  multiplied  by  a tri- 
gonometrical function,  Trauboth’s  fast  convolution  integration  algorithm  (11),  which  only  re- 
quiries  about  the  same  number  of  computation  steps  as  for  ordinary  integration,  is  used  here. 
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It  can  be  seen  from  Eqs.  (23)  and  (24)  that  each  successive  integration  advances  the  solution 
time  by  2r)  and  the  number  of  successive  integrations  required  depends  on  T and  n.  The  con- 
trol of  numerical  accuracy  lies  in  finding  the  zeros  of  and  the  subsequent  numerical  in- 
tegrations, and  both  are  well  established  numerical  techniques.  Other  quantities  such  as  (7*) 
can  be  found  by  the  same  procedure  or  calculated  from  their  relationships  with 

Equation  (20)  can  also  be  transformed  into  the  following  Volterra  integral  equation: 

r 

-n)  + I -2t,  -2t,  - r)dr 

- w°(.T  -r,)H{T  --n).  (25) 

On  close  examination,  however,  the  solution  of  Eq.  (25)  requires  repeating  exactly  the  same 
convolutions  as  in  Eq.  (24).  Iherefore  it  is  quite  unnecessary  to  use  the  numerical  techniques 
of  integral  equations,  even  if  only  one  single  integral  equation  is  involved,  for  obtaining  the  in- 
verse Laplace  transform  of  Eq.  (20).  On  the  contrary,  the  use  of  two  simultaneous  Volterra  in- 
tegral equations  for  obtaining  the  inverse  Laplace  transforms  of  and  w^,  as  in  Ref.  6 would 
have  unduly  increased  the  numerical  difficulty  and  the  computation  effort. 

Results  and  Discussions 

Numerical  results  are  obtained  for  a case  in  which  both  shells  are  made  of  steel  and  both 
fluids  are  water.  The  material  properties  and  dimensions  used  are  such  that 

f,  -p,  -i,  -C}  -17.79133 

M'  - 32.09875,  M - 6.41975 
h'la*  - 1/250,  h/a  - 1/50 

{ -0.8,  tj  -0.2.  (26) 

The  incidence  is  a step  wave  with  /(s)  — 1/s.  The  transient  responses  of  the  shells  are  calcu- 
lated for  four  transit  times  of  the  incident  wavefront,  i.e.,  T — 0 to  8.  For  this  duration  of  time, 
it  requires  twenty  successive  convolution  integrations  in  Eq.  (24).  The  parabolic  rule  is  used  for 
the  numerical  integration  in  conjunction  with  Trauboth’s  fast  convolution  scheme.  Since  nu- 
merical roundoff  and  truncation  errors  accumulate  at  each  successive  integration,  the  integra- 
tion intervals  At  are  kept  sufficiently  small  to  minimize  these  errors.  The  roots  of  A„  (s),  and 
from  these,  w^(T),  and  G„(r),are  evaluated  with  high  accuracy  before  starting  the  convolu- 
tion integration.  It  is  also  found  that,  similar  to  the  case  of  a single  shell  [7],  eight  terms 
(m  — 0 — 7)  in  the  series  solution  are  sufficient  for  the  representations  of  the  shell 
deflections,  velocities,  and  stresses.  The  integration  intervals  used  are:  At  — 0.01  for  wq  and 
w, , At  - 0.0025  for  »V2.  and  ” 0.00125  for  the  higher  terms.  If  the  solution  must  be  car- 
ried out  for  shorter  time  duration,  say  up  to  T — 4,  a much  coarser  At  will  suffice. 

In  all  subsequent  figures,  the  present  solutions  are  ploted  in  dotted  lines  and  compared  to 
the  solutions  obtained  by  the  method  of  Ref.  7 and  plotted  in  solid  lines  for  the  case  in  which 
the  outer  shell  is  absent. 

Figure  2 shows  the  results  for  the  wq's  and  wi’s.  Similar  to  the  single-shell  case,  the 
present  results  approach  their  respective  asymptotic  values  after  about  two  transit  times,  and 
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r-d-nw 

Fig.  2 - Time  histories  of  i»g  and  ui^ 

the  numerical  asymptotes  of  wq  and  •i’l  agree  with  formulae  (17)  and  (19)  respectively  within 
three  digits.  These  results  are  indicative  of  the  accuracy  of  the  present  solution  computation 
scheme.  For  this  particular  example,  the  outer  shell  is  rather  thin  and  its  effect  on  the  response 
of  the  inner  shell  is  observable  in  the  wg,  W],  and  wj  terms.  Its  effect  on  the  higher  modes  of 
the  inner  shell,  i.e.,  with  m > 2,  is  quite  undiscernible.  Sample  results  for  higher  w^’s  are 
exhibited  in  Fig.  3.  It  can  be  seen  there  that  the  present  results  for  wj  and  h>4  coalesce  with 
the  respective  results  for  the  outer  shell  absent  case.  The  same  can  be  laid  about  ^5,  and 
^7 , and  they  are  not  replotted  here. 
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Eight  terms  of  Wiji’s  arc  summed  in  Eq.  (7)  for  representing  the  total  response  of  the 
inner  shell.  Figure  4 plots  the  time  history  of  the  relative  radial  deflection  between  the  two 
apexes  of  the  inner  shell.  The  presence  of  the  outer  shell  causes  a downward  shift  of  the  curve 
representing  the  present  result.  This  is  due  to  the  diminished  wg  term. 

Figure  S shows  the  time  histories  of  the  radial  velocities  of  the  same  two  apexes.  The 
presence  of  the  outer  shell  causes  a slight  reduction  of  the  velocities  with  little  alteration  of 
their  profiles.  Again  this  can  also  be  expected  by  observing  the  slight  reduction  of  the  Wj  term. 


Fig.  4 - Relative  radial  deflection  between  the  Fig.  S - Time  histories  of  radial  velocities  at  the 
apexes  of  the  inner  shell  apexes  of  the  inner  shell 


It  is  reiterated  here  that  the  series  solution  in  the  form  of  Eq.  (7)  is  not  effective  for  ob- 
taining the  early  time  acceleration  and  pressure  if  the  incident  wave  has  a steep  front.  The  cus- 
tomary remedy  is  to  apply  a modified  Watson's  transform  to  Eqs.  (11)  through  (14)  and  calcu- 
late the  asymptotic  results  for  large  s [1].  This  is,  however,  much  more  complex  than  the  single 
shell  case  studied  in  Ref.  12,  and  is  outside  the  scope  of  the  present  report. 

The  polar  membrane  stress  (force  per  unit  length)  \f,  and  the  hoop  membrane  stress 
are  calculated  by 

■ (1  f.]Ml  +-»=»!  «l  <2’) 

and 

“Ti  f!)Mi  + -FMCOttll.  (28) 
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The  time  histories  of  the  polar  and/or  hoop  stress  ai  w — ir/2  and  at  « — »r  are  plotted  in 
Fig.  6 and  7 respectively.  The  shielding  of  the  outer  shell  reduces  the  stresses  slightly.  The 
mean  stresses  also  decrease  in  proportion  to  h’q.  The  maximum  stresses  occur  after  the  incident 
wave  has  completely  engulfed  the  outer  shell  and  are  about  twice  the  corresponding  static 


r - (T-Sr,)/{ 

Fig.  6 - Time  histories  of  the  polar  stress  at 
e = ir/2 


T - (T-nMi 

Fig.  7 - Time  histories  of  the  polar  and  hoop 
stresses  at  0 • tr 


All  results  of  the  present  example  infer  that  a thin  outer  shell  tends  to  be  transparent  to 
the  incident  pulse. 

In  conclusion,  it  can  be  said  that  an  accurate  solution  method  for  this  problem  has  been 
developed  and  can  be  used  for  parametric  studies.  Moreover,  a set  of  numerical  results  has 
been  meticulously  obtained  to  serve  as  a data  base  for  the  development  of  numerical  methods 
for  the  analysis  of  the  transient  response  of  double  hulls  of  practical  configurations.  One  im- 
mediate scheme  is  to  apply  the  currently  available  Doubly  Asymptotic  Technique  [2]  in  the 
outer  fluid  field  and  the  fluid  finite  element  technique  for  the  entrained  fluid. 
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NKF  ENGINEERING  ASSOaATES,  INC. 
8720  GEORGIA  AVENUE 
SILVER  SPRING.  MD.  20910 

ATTN  DR.  R.O.  BELSHEIM 

MAT  CO. 

2130  ARCH  STREET 
PHILADELPHIA,  PA.  19103 

ATTN  a McNAIGHT 


CATHOUC  UNIV.  OF  AMERICA 
DEa»T  OF  PHYSICS  AND  INSTITUTE 

OF  OCEAN  SCIENCE  AND  ENGINEERING 
WASHINGTON,  D.C. 

ATTN  H.  UBERALL 


CORNELL  UNIV. 

COLLEGE  OF  ENGINEERING 
DEPT  OF  THEORETICAL  AND 
APPLIED  MECHANICS 
ITHACA,  N.Y.  14860 

ATTN  Y.H.  PAO 


SANDIA  LABORATORIES 
SANDIA  CORPORATION 
DIVISION  9324 

ALBUQUERQUE,  NEW  MEXICO  87115 
ATTN  M.J.  FORRESTAL 
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